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Abstract 



^ . The use of the umbral formalism allows a significant simplification of the derivation of sum rules 

involving products of special functions and polynomials. We rederive in this way known sum rules 
and addition theorems for Bessel functions. Furthermore, we obtain a set of new closed form sum 
rules involving various special polynomials and Bessel functions. The examples we consider are 
relevant for applications ranging from plasma physics to quantum optics. 
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I. INTRODUCTION 



In this note we will show that the problem of deriving generating functions of special 
functions can be greatly simplified by the use of symbolic methods. This methods have been 
widely employed recently to deal with the Ramanujan master theorem l|-|3| and to derive 
closed expressions for integrals involving various combinations of Bessel functions. 

The methods we are going to describe in the paper complete previous researches summa- 
rized in Ref. where it was shown that the use of operational techniques allowed us the 
derivation of sum rules hardly achievable with conventional means. 

Before entering into the specific details of the paper we will premise the main ingredients 
of the formalism we will employ in the following, to establish known and new sum rules 
within simplified and unifying point of view. We start by introducing the umbral notation 
with the operator c: 

C"^0 = ifa, = T, (1-1) 

— c^), without any restriction on a. 



According to the previous definition, Bessel functions of the first kind jsj] can be written 



as 



2k+a 



k=0 



(-l)Mx/2) 



k\T(a + k + l] 



C-) e 



-cix/2f 



:i.2) 



while the Tricomi-Bessel functions of order a [5|-[7| are given by 

\k 



CJx) 



fc=0 



-XI 



k\T(a + k + l) 



0£i 

= x-"/V«(2v^) = c"e" 



:i.3) 



For the two-variable Laguerre polynomials 



Big 



9|, the same notation yields. 



Ln{x,y) =n\Y^ 



k„,n—k 



fc=0 



X 



:i.4) 



where Ln{z) is the ordinary Laguerre polynomial. Finally, for the Wright function 10| we 
have 

^-(^1/^) = E ITFTT^-TX = ^"'^ e'^'^^o. (1.5) 



^ r! Tiu + fir) 
Given a function with series expansion 



X 



:i.6) 



n=0 



we will define the corresponding Laguerre-based function as 

oo 

Lf{x,y) = y2^^n{x,y). (1.7) 

n=0 

In an analogous way, we define the Hermite-based functions as (m G 1^+) 

oo 
n=0 

with 

[n/m] 

{n — mk)\ k\ 



Ht>(x,y)=„'.j:j—J-^. (1.9) 



fc=0 

For future reference, we remind that the two-variables Hermite polynomials verify the iden- 
tity 

Hjr\-x,y) = {-irHjr^{x,{-iry) (I.IO) 
and their generating function is given by 

Y.-Hjr\x,y) = e^^^y'\ (1.11) 



n=0 



The results we will obtain in this paper will also be exploited to derive families of sum rules 
useful in classical and quantum optics. In particular, we will make a comparison with the 



conclusions of the recent analysis of Ref. 



Ill ], in which the sum rules for the product of 



Bessel functions have been used to derive the expression for the work done by a modulated 
electric field forcing a damped harmonic oscillator. 

II. GENERATING FUNCTION OF BESSEL FUNCTIONS AND SPECIAL POLY- 
NOMIALS 

We prove the efficiency of the procedure by means of some introductory examples, which 
are a direct consequence of the previous definitions. We start with a well-known result (see 



Eq. (5.7.6.1) from 12|). By using Eq. (11. 2p . it is a simple task to show that 

t"" ( x\^ ( c 1 

2^:^/n+u{x) = i^c-j exp j -- (x^ - 2xt) ^ V?o 



n=0 

u/2 



^^^^j MVx^-2xt), {\2t\<x, z/gC), (2.1) 



obtained here in a nonconventional and very straightforward way. 
Furthermore, one has 



n=0 ■ ^ ' 



where 



k=0 



(2.2) 



(2.3) 



are the so-called Hermite-based Tricomi functions 



9|. To get the above sum rule, Eqs. fll.lOp 



and (11. lip have been used, along with the following identity (m G 



n=0 



^0 = nc^r^ {x,{-iry) 



(2.4) 



that proves to be useful also to show that 



E^^"M(^)=^^o"'K2 

n=0 ■ ^ 

where we have introduced the new functions 



(m)A/a;\V'" X 



m 



(2.5) 



°° tt('"t-) / \ 



k=0 



k\T{iik + iy + l] 



(2.6) 



i.e., the Hermite-based functions obtained for a„ = l/r{fin + u + 1). 

By using the result (ll.4p we can easily derive the new (hardly achievable with conventional 
means) sum rule 



^ — Jn{z) Ln{x, y) = lCqI- 
n=0 ^' ^ 



xt z z{z — 2yt) 



(2.7) 



where the Laguerre-based function is obtained from Eq. (12. 3p by replacing Hn"^^ with 
Ln- In the same way we can obtain a generating function involving the product of Hermite 
and Laguerre polynomials, which can be written as 

J] - L„ (x, y) Hj^^ iz,w) = ey' ^^+y hC^^ {xt{z + 2ywt),x^ wt") , (2.8) 



n=0 



through the use of Eq. (12. 4p . 



III. SYMBOLIC DERIVATION OF SOME ADDITION THEOREMS 



Our method allows us to prove the Graf's addition theorem [5| in a very straightforward 



a 



way. By using Eq. (11 ■2p and the identity 12 1 

Uz) = exp U (t-l)], (3.1) 

n=— oo ^ ^ ^ 

we get the formula {u E C) 

oo 

G^''\x,y;t) = J2 t^JnUx)Jn{y) = (c|)' e-^'-^'-\o (3.2) 



where 



t = - [x —xyt) Y = T — • 

^ 4 ^ ^ ^ ^ Axt 



After a simple manipulation, we can write 



G(^^(a^,y;t)= (^)'f:^^|^^.-.(2v^). (3.3) 

and, since (see Appendix) 

V^^J,_„(x)=P^ MVx^-2xt), (|2t|<x, z/eC), (3.4) 

n! \ X J 



n=0 

one obtains 



g^'^H3:,y;^)= ( "'^J'^^ J ' J.(Vx^+y^-Xy(t + t-l)), (a;>^/t). (3.5) 

Then, setting t = e*^, the Graf's addition theorem results: 

J2 e'^' Jn+u{x) Jn{y) = (^^^^^] J.(A/x2 + y2_2x2/cos^). (3.6) 

n=-oo ^ ^ ^ ^ 

In the same way, we can evaluate the following extension of the Neumann addition the- 
orem 

E rj.(x) J,.(y) = (x' ^11^ ) ' ^3.7) 
where we have introduced the hybrid special functions defined with (12. 3p as: 

^i^(™)(x,y|0 = ^Cil+,(x,y) (/i e M). (3.8) 

A;=0 



In spite of their awkward aspect, these functions have quite simple properties and their 
study may shine further hght on the properties of multivariable Bessel functions, which 
are playing a role of central importance in the theory of multi-photon processes induced 



19| . The same procedure can be applied to the 



by charged systems in intense laser fields 
derivation of sum rules involving products of three Bessel functions. 
Also the evaluation of sums of the type 111 



(3.9) 



required in some problems in plasma physics and quantum optics, can be performed quite 
straightforwardly using our symbolic method. By noting that 



(3.10) 



i?=0 



from the Graf's addition theorem one has 



J2 e^'"'' Jn+i{x) Uy) 



x — ye^^\ ( ^x^ + y"^ — 2xy cos'd 



exp < — c 



y^o, (3.11) 



and by applying to this expression the so-called Hoppe formula for the derivative of composite 



functions 



with 



id 



dt 



fc=0 



k 



k\ 



(3.12) 



^rn,.(/W) = E 



Jm 



we finally obtain 



j=Q V"' / k=0 ^ ^ p=0 ^' 



(3.13) 



where 



5=0 



P 



1 \ / q 



2 I ^ \r 

r=0 



E P-^) 



X 



k+p ( _yy-k+p 



{x - yY+P 



Ji+p{x-y). (3.14) 



In a similar manner, by using Eq. (12. ip specialized to the case of integer z/, we can write 
(m, / G Z) 



n=0 



X] 



. a; 
c — 
2 



exp 



n=0 

x"^ — 2xt 



i=l 



(3.15) 



i=l 



and, by recalling that 14l-ll6| 

m 

{tdtr = Y,S2{m,k)t^d[ 

k=l 

with S2{k,'m) Stirling numbers of the second kind, we find (see also Ref. [l 
Et\x) = J2S2{m,k) (c|) 

k=l 



(3.16) 



fc=l 



X — 2xt 



(3.17) 



where C/+fc is the Tricomi-Bessel function defined in Eq. (11. 3p . 
All the new results obtained above have been verified numerically. 

In this paper we have shown that using a symbolic procedure of umbral type, significant 
progress can be achieved in the theory of sum rules involving Bessel functions. Such a point 
of view has opened a new scenario offering a wealth of analytical possibilities amenable for 
further generalizations of practical interest in various fields involving the theory of special 
functions. We have indeed mentioned laser assisted processes in strong laser fields 19|, but 
we should also include synchrotron radiation problems regarding the emission of electron 
radiation in nonconventional undulator magnets 20|. Within this last framework a signif- 



icant role has been played by the so called multi-index Bessel functions 
can be greatly simplified by the umbral formalism discussed here. 



2l|, whose theory 



IV. APPENDIX 



In this Appendix we prove the identity (13. 4p . used in sec. 
theorem. By starting from the familiar derivative formula |5| 



TT]to prove the Graf's addition 



{x-'d,y [X" Ux)] = X"-^ J,.n{x) 



[u G C,n = 0,1,2,...), 



(4.1) 



we obtain 

exp{TX-'d.,} [x'^J^{x)]=x''J2^^J^-k(^) (4-2) 



and, by using the operational identity 



22| 



fc=0 



exp {tx-' 4} fix) = /(Vx2 + 2r), (4.3) 

we get 

f (l^ . ( :/EIliy .uVx^rj-r) (4.4) 

k=0 ' \ ^ J 

from which, after setting r = —xt, we recover (13.41) . 
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